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ABSTRACT 

Interactions with neutrons and protons play a crucial role for the neutrino opacity of matter in the supernova 
core. Their current implementation in many simulation codes, however, is rather schematic and ignores not 
only modifications for the correlated nuclear medium of the nascent neutron star, but also free-space correc¬ 
tions from nucleon recoil, weak magnetism or strange quarks, which can easily add up to changes of several 
10% for neutrino energies in the spectral peak. In the Garching supernova simulations with the Prometheus- 
Vertex code, such sophistications have been included for a long time except for the strange-quark contributions 
to the nucleon spin, which affect neutral-current neutrino scattering. We demonstrate on the basis of a 20 M© 
progenitor star that a moderate strangeness-dependent contribution of gl = -0.2 to the axial-vector coupling 
constant ga ~ 1 -26 can turn an unsuccessful three-dimensional (3D) model into a successful explosion. Such a 
modification is in the direction of current experimental results and reduces the neutral-current scattering opac¬ 
ity of neutrons, which dominate in the medium around and above the neutrino sphere. This leads to increased 
luminosities and mean energies of all neutrino species and strengthens the neutrino-energy deposition in the 
heating layer. Higher nonradial kinetic energy in the gain layer signals enhanced buoyancy activity that en¬ 
ables the onset of the explosion at ~300ms after bounce, in contrast to the model with vanishing strangeness 
contributions to neutrino-nucleon scattering. Our results demonstrate the close proximity to explosion of the 
previously published, unsuccessful 3D models of the Garching group. 

Subject headings: supemovae: general — hydrodynamics — instabilities — neutrinos 


1. INTRODUCTION 

According to the standard paradigm of the explosion mech¬ 
anism of core-collapse supernovae (SNe), neutrino heating 
above the gain radius initiates the outward acceleration of 
the stalled bounce shock and provides a major fraction of 
the energy that unbi nds the explosion ejecta and powers the 
SN blast wave (e.g . , [Colgate & White|1966 JBethe & Wilson 

IJankall 


19851 [Bethe||1990 |janka||2012t Burrows||2()13| ). Except for 

stars near the low-mass end of SN progenitors with oxygen- 
neon-magnesium cores or small iron cores and a steep den¬ 
sity gradient of the overlying shells, successful explosions 
cannot be obtained in spherical symmetry (ID) by state-of- 
the-art simulations. Eor more massive progenitors the suc¬ 
cess of the neutrino-driven mechanism requires the support by 
multi-dimensional hydrodynamic flows in the postshock layer 
associated with convective overturn and the sta nding accre¬ 
tion shock instability (SASI; [BTondin et al.|20Q3| ). Such flows 
increase the neutrino-energy deposition and creke buoyancy 
and turbulent pressure, thus reducing the critical threshold for 
the neutrino lu minosity to trigger the explosive runaway of 


& Muller 19961 [Murphy & Burrows[2008 

[Nordhaus et al. 

2010 Hanke et al.[[20I2[ Dolence et al.[[201 

3[ [Murphy et al. 


taneous mass accretion and outflow from the proto-neutron 
star (PNS) after shock revival, which steepens the rise of the 
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explosion energy ( Marek & Janka|200^ . 

The first self-consistent three-dimensional (3D) stellar 
core-collapse simulations with successful explosions were 
performed by |Eryer & Warren (2002 2004| ), using smoothed- 
particle hydrodynamics and a simple, grey diffusion descrip¬ 
tion of neutrino transport, which favored rapid explosions 
in 2D and 3D. Recently, the Garching group has obtained 
a 3D explosion for a 9.6 M© star with sophisticated, state- 
of-the-art ray-by-ray-plus (RbR-f), multi-group transport and 
found that 3D hydrodynamic flows enhance the explosion en 


ergy c ompared to the 2D (axisymmetric) case ([Melson et al. 


2QI5| ). However, previous 3D simulations of this group for 


11.2, 20, and 27 M© progenitors (covering -400-550 ms af¬ 
ter bounce) did not produ ce explosions although correspond¬ 
ing 2D models exploded ( |Hanke et al.||2013[ |Tamborra et al. 


[MT3l |2014a|b| ). This is consiste nt with studies based on 


parametrized neutrino heating (e.g . ,|Hanke et al.|2012 


Couch 


|2013t [Couch & Q’Connor||2014t [Couch & Ott||20i5| ) and 
self- consistent simulations with approximate neutrino trans¬ 
port ( [Takiwaki et al.[20T2 2014 ), which revealed that the on¬ 
set of explosion is less faWed or delayed in 3D com pared 
to 2D, contradicting claims by Nordhaus et al.[ (|2010[) and, 
though moder ated, by Burrows et al.[ ( |2012[ ) and [Dolence 


et al. 


([2013|). Also [Lentz et al.[ ([2015|), using high-fidelity RbR 
multi-group neutrino diffusion, highlight a successful 15M© 
explosion that sets in significantly later than the correspond¬ 
ing 2D model. Large-scale deformation modes seem to have 
the strongest supportive effect on the explosion but are weak¬ 
ened by the forward turbulent energy cascade to small-scale 
structures in 3D, which is opposite to the 2D case. 

While the resolution feasible in current full-scale 3D SN 
models is st ill insufficient to satisfactorily represent postsh ock 
turbulence ( [Abdikamalov et al.|20I4t[Radice et al.|20I5[ ), the 
results reported above suggest that important physics might 
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Fig. 1.— 3D iso-entropy surfaces for different times after bounce for exploding model 3Ds. Colors represent radial velocities. The supernova shock is visible 
by a thin surrounding line, the proto-neutron star by a whitish iso-density surface of 10^^ gcm“^. The yardstick indicates the length scale. Strong SASI activity 
occurs between ~120ms and ~280ms. (An animation and interactive version of this figure are available in the HTML version of this article.) 


Still be missing in the models. One of the aspects to be scruti 
nized are the pre-collapse initial conditions, which result from 
ID stellar ev oluti on modeling. |Couch & Ott| ( |2Q13| ), |Couch 


|et al.| ( |20T5] ), and |Muller & Janka| ( |2015| ) indeed confirmed 
speculations that large-amplitude perturbations of low-order 


modes in the convective shell-burning layers (e.g., | Arnett & 
|Meakin|2011 and references therein) might facilitate the de¬ 
velopment of explosions. Further progress will require 3D 
modeling of the final stages of stellar evolution. 

Here we demonstrate that remaining uncertainties in the 
neutrino opacities, in particular the neutrino-nucleon interac¬ 
tions at subnuclear densities, can change the outcome of 3D 
core-collapse simulations. As an example we consider pos¬ 
sible strange-quark contributions to the nucleon spin in their 
effect on weak neutral-current scatterings. We show that a 
moderate isoscalar strange-quark contribution of gl = -0.2 to 


the axial-vector coupling constant ga = 1-26, which is not far 
from current experimental results, suffices to convert our pre¬ 
vious non-exploding 20 M© 3D core-collapse run into a suc¬ 
cessful explosion. 

We briefly describe our numerical approach in Sect.j^ sum¬ 
marize basic facts about the strange-quark effects in neutrino- 
nucleon scattering in Sect. discuss our results in Sect. 
and conclude in Sect. [5j 

2. NUMERICAL SETUP AND PROGENITOR MODEL 


We performed 2D and full (47r) 3D simulations of a nonro- 
tating, solar-metallicity 20 M© pre-SN progenitor (|Woosley & 
|Heger|2007| ). 

We used the Prometheus-Vertex hydrodynamics code with 
three-flavor, energy-dependent, ray-by-ray-plus (RbR-f) neu¬ 
trino transport including the full set of neutrino reactions and 
microphysics ( Rampp & Janka||2002[ [Buras et al.||2006 ) ap- 
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Fig. 2.— Cross-sectional entropy distributions (in per nucleon) for the 3D models without (3Dn; upper row) and with strangeness contributions (3Ds; 
bottom). The bottom row clearly shows stronger SASI activity in model 3Ds (180 ms, 250 ms), whose traces are still imprinted on the ejecta geometry after the 
onset of the explosion (530 ms; note the different scale). 


plied in 3D also by |Hanke et~ar ( [2Q13| ) and |Tamborra et al.| 
( 2013 |2014a|b|), in particular the high-d ensity equation of 
state (EoS) oi |Lattinier & Swesty| ( |1991| ) with a nuclear in- 
compressibility of K = 220 MeV. To avoid time-step con¬ 
straints, the convectively stable central core with a radius of 
10 km was computed in spherical symmetry. We used a ID 
gravity potential (which is unproblematic for strongly PNS- 
dominat ed gravity fields) in cluding general relativistic cor¬ 
rections ( |Marek et al.|2006| ). The radial grid had a reflecting 
boundary condition at the origin and an inflow condition at the 
outer boundary of 10^ cm. It had 400 nonequidistant zones 
initially and was reflned in steps up to >500 zones, provid¬ 
ing an increasingly better resolution of Ar/r ~ 0.01... 0.004 
around the gain radius. For the neutrino transport 12 geomet¬ 
rically spaced energy bins with an upper bound of 380 MeV 
were used. The growth of non-radial hydrodynamic instabili¬ 
ties was seeded by imposing random cell-to-cell perturbations 
of 0.1% in density on the whole computational grid 10 ms af¬ 
ter bounce. 


3. STRANGENESS CONTRIBUTIONS TO NEUTRINO-NUCLEON 
SCATTERING 

The lowest-order differential neutrino-nucleon scattering 
cross section reads 

^ = ^^[c^(l+cos0) + c^(3-cos0)] , (1) 

with 6 being the incoming neutrino energy, 6 the scattering 
angle, Gp Fermi’s constant, and Cy and vector and axial- 
vector coupling constants, respectively. The latter are Cy = 
\ - 2 sin^ ^ 0.035, 0.63 for yp vp and 

Cy = = “^a/2 ~ -0.63 for vn vn with ^ 1.26 

and sin^^w ~ 0.2325. For iso-energetic scattering {e' = e), 
Eq. ^ yields the total transport cross section 

. C dCTn 2Gp6^ / ^ 


While in our SN simulations corrections due to nucleon ther¬ 
mal motions and recoil, weak magnetism, and n ucleon cor¬ 
relations at high densities are taken i nto account ( Rampp &| 
Janka|200'2{ |Buras et al.|2006] ), Eqs. ( |1|2| ) provide good esti¬ 


mates. Strange quark contributions to the nucleon spin modify 
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Fig. 3.— Explosion diagnostics for model 3Ds (thick lines) compared to the non-exploding model 3Dn (thin lines) as functions of post-bounce time ^pb. Top 
left: Angle-averaged shock radius (black), gain radius (red) and NS radius (blue; defined by a density of 10^^ gcm“^); top right: diagnostic energy (positive 
total energy behind the shock). Gray lines display the corresponding 2D models without (2Dn, thin) and with strangeness contributions (2Ds, thick); middle 
left: mass-accretion rate (M) ahead of the shock (red) and baryonic NS mass (blue); middle right, bottom left and right: mass, non-radial kinetic energy, and 
time-integrated neutrino-energy deposition in the gain layer, respectively. 


Ca according to 


= 2 (±ga 




( 3 ) 


where the plu s sign is for vp and the minus sign for vn scatter¬ 


ing (s ee, e.g., [Horowitz 


Langanke & Martmez-Pinedo 


2003]). Since < 0, the cross section tor y/7-scattering is 
increased and for y^-scattering decreased. 


Employing Eq. ([2|) with gl = -0.2, Horowitz (2002) es¬ 
timates 15, 21, 23% reduction of the neutral-current opac¬ 
ity for a neutron-proton mixture with electron fractions = 
0.2, 0.1, 0.05, which are typical values for the layer be¬ 
tween neutrinosphere (at density p ~ 10^^ gcm“^) and p ~ 
10^^ g cm“^ for hundreds of milliseconds after bounce. Since 
strangeness does not affect charged-current interactions and 
NS matter is neutron-rich, the reduced scattering opacity al¬ 
lows mainly heavy-lepton neutrinos (Vx = v^, v^, Vr) to 


leave the hot accretion mantle of the PNS more easily. This 
was found to enhance the expansion of the stalled SN shock 
in ID models, although not enough for successful shock re¬ 
vival ( jETebendorfer et al.|2002] [Langanke & Martmez-Pinedo 
[200^ . However, below we will show that the situation can be 

fundamentally different in 3D simulations. 


4. RESULTS 

We compare 2D and 3D core-collapse simulations of the 
20 Mq star with strangeness corrections in neutrino-nucleon 
scatterings, using gl = -0.2 (models 2Ds, 3Ds), to corre¬ 
sponding simulations without strange quark effects (g| = 0; 
models 2Dn, 3Dn) as in all SN simulations of the Garching 
group so far. To explore “extreme” effects, our choice of gl 
is by its absolute value somewhat bigger th an theoretical and 
experimental determinations of gl 0.1 ([Ellis & Karliner 


1997[[Alexakhin et al. [2007[[Airapetian et ar.|2007[). 
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Neutrino-driven explosion of 20 solar mass progenitor in three dimensions 
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Fig. 4.— Neutrino-related quantities for model 3Ds (thick lines) compared to the non-exploding model 3Dn (thin lines) as functions of post-bounce time ^pb- 
Top and middle left: Neutrino luminosities and mean energies (defined as ratios of neutrino energy and number fluxes), respectively, for Ve (blue), Ve (red), and 
one kind of Vx (black) for a distant observer at rest; bottom left: volume-integrated neutrino-energy deposition rate in the gain layer (upper panel) and cooling 
rate between NS radius and gain radius (lower panel); top, middle, bottom right: advection (blue) and heating time scale (red), corresponding time-scale ratio, 
and heating efficiency in the gain layer, respectively. 


4.1. Dynamical evolution towards explosion 

The sequence of 3D images in Fig.shows the post-bounce 
evolution of the exploding model 3Ds; the entropy cuts in 
Fig. [^demonstrate the differences to the unsuccessful model 
3Dn. In both cases the dynamics of the accretion layer are 
strongly SASI-dominated. In model 3Ds there is no indica¬ 
tion of postshock convection before SASI sloshing motions 
first appear at post-bounce (p.b.) time ^pb ~ 120 ms. These 
reach full strength around 180 ms p.b. and continue in varying 
directions until ~280ms. Only later on convective overturn 
takes over as the dominant non-radial instability. In model 
3Dn moderate buoyancy activity is visible from 100-180 ms 
before SASI becomes dominant, too. Early convection is en¬ 
abled by a slightly larger radius of the stalled shock in 3Dn 
as a consequence of a slightly lower mass-accretion rate M 
(Fig.§. This difference improves the growth conditions for 


postshock convection. In 3Dn an erroneous change of the 
transition between low-density and high-density EoS caused a 
delay of the core collapse. Therefore the mass-accretion rate 
until r ~ 150 ms is slightly reduced and the Si/Si-fO interface 
arrives at the shock ~15 ms later. These early differences are 
inessential for our discussion because neutrino heating creates 
favorable conditions for the explosion of model 3Ds only later 
than -300 ms. 

From ^pb ~ 170 ms on, 3Ds exhibits clearly larger SASI 
amplitudes and higher postshock entropies, which increases 
the maximum and average shock radii (Figs.[^[^. This model 
also shows larger non-radial kinetic energies m the gain layer. 



(Fig. [ 5 , except during ^pb - 200-300 ms, when a powerful 
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Fig. 5.— Total mass-accretion rate in inflows (v^ < 0; bottom) and mass- 
outflow rate in rising plumes in the gain layer (v^ > 0; top) at three times after 
the onset of explosion of model 3Ds. Vertical dashed lines mark the positions 
of the average gain radius, the steep forward edges indicate the average shock 
radii. 

spiral SASI mode develops in model 3Dn but not in 3Ds, al¬ 
beit without pushing the shock s ufficiently far out for revival 
(in conflict with recent results by |Fernandez|2015 ). 

At fpb > 300 ms, roughly 50 ms after the Si/Si-fO inter¬ 
face has fallen through the shock and the shock has expanded 
to -150 km, the evolutions of models 3Dn and 3Ds separate. 
While in 3Dn the average shock radius, (/?shock), retreats again 
and conditions become unfavorable for an explosion, model 
3Ds exhibits positive trends in all explosion-diagnostic pa¬ 
rameters (e.g., (/?shock), ^gain)- Continuous shock ex¬ 

pansion signals runaway and finally outward acceleration sets 
in at fpb > 360 ms, at which time the recombination of free 
nucleons to a-particles in the largest plumes begins to release 
energy and the instantaneous “diagnostic energy”, 

^exp = f dV p^tot, (5) 

»^etot>0,postshock 

starts to rise steeply, correlated with a fast growth of the mass 
in the gain layer (Fig.[^. In Eq. ^ the volume integration is 
performed over the postshock region where the total specific 
energy, 

etot = e+ I|v|^ + O + [ebind(^®Fe) - Cbind] , (6) 

is positive, with e, ^|vp, and O being the specific internal, 
kinetic, and (Newtonian) gravitational energies. The brack¬ 
eted term expresses the difference between the specific nu¬ 
clear binding energy when all nucleons are finally recombined 
to iron-group nuclei compared to the nuclear composition at a 
given time. It therefore accounts for the maximum release of 
nuclear energy and corresponds to an upper limit of £’exp. 

Both corresponding 2D models, 2Dn and 2Ds, also explode 


after the Si/Si-fO interface has passed the shock, but in 2Ds 
the outward shock acceleration and rise of £’exp sets in -50 ms 
earlier (Fig.[^. Strangeness corrections therefore create more 
favorable explosion conditions also in the 2D case, although 
their influence is modest in successful models, similar to their 
small effect in the ID cas e, which is far away from explosion 
( [Liebendorfer et al.|2002| ). 

4.2. Strangeness corrections and explosion 

The impact of strangeness effects on the neutrino emission 
and the explosion is displayed in Fig. Model 3Ds con¬ 
sistently exhibits higher luminosities and mean energies of 
the emission for all neutrino species, and, consequently, a 
higher neutrino-energy deposition rate in the gain layer, Qgain, 
a higher heating efficiency. 


gain 2gain 


(7) 


a smaller gain radius (Fig.[^, and a shorter heating time scale. 


Theat 



( 8 ) 


with ^eain = &Vp{e + l\v\^ + <I)) being the binding 

o J\/<gain/ ^ 

energy of the gain layer. Since the effective time scale of mass 
advection through the gain layer. 


Tadv 


M, 


gam 


M 


(9) 


(where M > 0), which measures the average exposure time 
of matter to neutrino heating, is very similar in models 3Ds 
and 3Dn, the smaller Theat in 3Ds also leads to a higher time- 
scale ratio Tadv/Theat- The ratio Tadv/Theat exceeds the critical 
value of unity shortly before the SN shock in 3Ds begins its 
runaway expansion. 

The mean energies of the radiated neutrinos in model 3Ds 
are up to -1 MeV higher and the luminosities of and 
by up to -10-15%, whereas the Vx-luminosities rise by up to 
-30%. The increase of the total neutrino luminosity is more 
than 6 x 10^^ erg s“^ at maximum, which mainly comes from 
layers below the v^-sphere, because the neutrino-loss rate 
2cooi between the location of this sphere (at -lO^^gcm"^) 
and the gain radius differs between models 3Ds and 3Dn by at 
most -10^^ erg s“^ (Fig.[^. Note that at ^pb > 300 ms the rel¬ 
ative differences of the neutrino properties of models 3Ds and 
3Dn decrease and even change sign, because the former ex¬ 
plodes whereas the latter continues to collapse and to accrete 
mass onto the PNS at a higher rate. 

For Ye = 0.1-0.05, strangeness effects in the neutrino- 
nucleon interactions reduce the effective opacity /Ceff = 
V^abs(kabs + kscatt) for Only by 2-3% and for by 8-10%. 
A considerable part of the observed luminosity enhancement 
of Ve and is therefore caused indirectly by a stronger con¬ 
traction of the PNS in response to the larger energy loss 
through Vjc emission. With the smaller PNS radius (Fig. 
and steeper density profile, the neutrinospheres of and Vg 
move inward to higher temperatures. 

The reduction of the weak neutral-current scattering by 
strange quark contributions to the nucleon spin therefore en¬ 
hances the neutrino luminosities and mean energies directly 
and indirectly and thus strengthens the neutrino heating in the 
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gain layer. This amplifies buoyancy and turbulent mass mo¬ 
tions behind the shock, which is signaled by higher non-radial 
kinetic energy in model 3Ds at ^pb > 300 ms (Fig.l^, thus fos¬ 
tering the explosion of this model in contrast to 3Dn. 


4.3. Further development of explosion energy 

At the end of the simulation at 530ms p.b., the PNS has 
a baryonic mass of 1.91 M© and the shock has expanded to 
an average radius of (T^shock) ~ 1000 km. The diagnostic en¬ 
ergy, £’exp, has reached 0.2 x 10^^ erg and rises steeply and 
linearly with a rate of ^exp ~ 1.2 x 10^^ erg s“^. This growth 
rate can be understood by the ejection of a continuous outflow 
of freshly neutrino-heated matter, Mv,>o ~ 0.08-0.1 M© s“^ 
(Fig.[^ the subscript indicates positive radial velocity), which 
is mainly fed by the shock-accreted matter that is channeled 
to the gain radius in persistent accretion downdrafts (Fig. 
with mass-inflow rates up to -0.3 M© s“^ (Fig.|^. The mass 
outfiow does not only absorb neutrino energy, it also releases 
nuclear binding energy from the recombination of neutrons 
and protons to a-particles and heavy nuclei. Since neutrino 
heating rough ly neutralizes the gravitational bin ding energy 
of the matter ( |Janka|20QT| Marek & Janka|200^ , the growth 


of £’exp can be estimated in terms of Mv^>o and an average 
nuclear recombination energy per nucleon, 6nuc, as 


^exp 


M. 


^nuc 

■>0 - 

me 


1.2 X 10^* 


/ ^Vr>G \ / ^nuc \ 

1,0.1 Mo/s/l6MeV/ 


( 10 ) 

To unbind the overlying layers of the progenitor with a grav¬ 


itational binding energy of -0.6x10^^ erg would therefore re¬ 
quire massive outfiow for several 100 ms. The neutrino-driven 
wind after the end of accretion as well as explosive nuclear 
burning will provide additional energy. For a reliable deter¬ 
mination of the final SN explosion energy, model 3Ds would 
need to be evolved considerably longer than in our simulation. 


5. CONCLUSIONS 

We showed that strangeness contributions to neutrino- 
nucleon scattering with an axial-vector coupling of gl = -0.2 
are sufficient to turn a non-exploding 3D simulation of a 
20 Mq model (in which ga = 1 .26 was used for the standard 
isovector form factor) to a successful explosion. Strange- 
quark effects in the nucleon spin reduce the neutrino opac¬ 
ity of neutron-rich matter inside the neutrinosphere and thus 
directly and indirectly enhance the luminosities and mean 
energies of the radiated neutrinos. This leads to amplifi¬ 
cation of the neutrino-energy deposition behind the stalled 
shock because charged-current interactions are not affected 
by nucleon strangeness. Owing to a reduced neutrino-heating 
time scale and stronger non-radial mass motions in the gain 


layer, the shock is driven to runaway expansion -100 ms af¬ 
ter the passage of the Si/Si-rO interface. The enhanced neu¬ 
trino emission enabled by the strangeness effects is associ¬ 
ated with (and partly caused by) a faster contraction of the 
PNS and a corresponding rise of the neutrino spheric tem¬ 
peratures. Strangeness corrections in neutrino-nucleon scat¬ 
tering are therefore similarly beneficial for shock revival as 
“softer” nuclear EoSs, which also lead to faster PNS con¬ 
traction and stronger emission of more energetic neutrinos 
dMarek & Janka|2009[|JMika|2012[[Suw et al.|2013|l. 


Our results demonstrate how close previous, unsuccessful 
3D core-collapse models with the Prometheus-Vertex code 
were to explosion. They also underline that an accurate 
knowledge of neutrino-nucleon interaction rates, in particu¬ 
lar also for neutral-current scattering, is of crucial importance 
for assessing the viability of the neutrino-driven explosion 
mechanism. While strangeness contributions affect neutrino- 
nucleon scattering everywhere, such opacity modifications 
between the subnuclear regime and neutrinospheric densities 
are most relevant for SN shock revival. In the discussed 20 M© 
simulations a modest -15% diminution of neutrino-neutron 
scattering makes all the difference between explosion and fail¬ 
ure. Theoretical and experimental determinations yield some- 
what smaller absolute values for gl than assumed in our study 


( Ellis & Karliner|1977}|Alexakhin et alJ2007t [Airapetian et^ 
|2007| ). However, in-niedium effects like coirelations in low- 
density nucleon matter may cause similar opacity reductions 
(C. Horowitz, private communication). 

Nucleon-strangeness effects should also be investigated in 
3D simulations of other progenitors. Moreover, our calcu¬ 
lations must be repeated with better than 2° angular zoning 
to ensure that the explosion is robust and withstands higher 
angular resolution of the cascading of turbulent energy from 
the largest scales to small structures , which can be harmful 
for shock revival ([Hanke et al.||2012||Couch||2013||Couch~^ 


Q’Connor|2Q14[[Abdikamalov k al.|2Q14| ). In any case, how- 

ever, the outcome of multi-dimensional core-collapse simu¬ 
lations that marginally explode or fail can sensitively depend 
on effects on the 10% level in the neutral-current neutrino- 
nucleon interactions. This sensitivity needs to be taken into 
account in numerical implementations of these rates and 
might also be important for understanding partially conflict¬ 
ing model results published by different groups. 
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